Non-Abelian geometrical control of a qubit in an NV center in diamond 
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We propose an approach for an optical qubit rotation in the negatively charged nitrogen-vacancy 
(NV) center in diamond. The qubit is encoded in the ground degenerate states at the relatively low 
temperature limit. The basic idea of the rotation procedure is the non-Abelian geometric phase in 
an adiabatic passage, which is produced by the nonadiabatic transition between the two degenerate 
dark states. The feasibility is based on the success of modeling the NV center as an excited-doublet 
four-level atom. 
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I. INTRODUCTION 



Phases play a major role in all interference and diffrac- 
tion phenomena in optics, wave physics and quantum me- 
chanics. The study of phase has generated a large num- 
ber of practical device [l|: superconducting interference 
device, interferometers to measure lengths and small ro- 
tations, modulators to modulate wavefront, etc. In the- 
oretical physics, phases are involved in the distinction 
between fermions and bosons, decoherence of a quantum 
system and Aharonov-Bohm effect. It is well known that 
a quantum system acquire a Berry phase Q which evolves 
adiabatically around a circuit by varying parameters in 
its Hamiltonian. This phase is a sum of two parts. The 
first one is related to the instantaneous energy of the 
system, which is called dynamical phase. The second 
one is a purely geometric property of the closed path fol- 
lowed by the system in space, which is called geometrical 
phase. Later, non-Abelian gauge fields arising in the adi- 
abatic development of simple quantum mechanical sys- 
tems are introduced by Wilczek et al [3j. Nowadays, the 
Berry phase is not only of interest from a fundamental 
point of view, for example, the induced gauge field [1, [B[ 
and topological states of matter Q, but also may have 
important applications in quantum information process- 
ing. Holonomic quantum computation Q is a prototype 
of research on quantum gates based on Abelian or non- 
Abelian geometric phases. Its significant attraction lies 
on that these quantum gates could be inherently robust 
against some local perturbations. Moreover, interest are 
paid on the adiabatic evolution of dark eigenstates of the 
considered system in order to remove any accompanied 
dynamical phase shift. 

On the other hand, the negatively charged nitrogen- 
vacancy (NV) center in diamond has recently emerged 
as a promising candidate for practical and scalable im- 
plementation of quantum information processing due to 
the following robustness: it is an individually addressable 
quantum system; its quantum state can be initialized 
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and manipulated with optical and microwave fields, and 
even measured with high fidelity at room temperaturefH- 
[l2j . NV centers are formed by a substitutional nitro- 
gen atom adjoining to a vacancy in the diamond lattice. 
The dangling bonds near the vacancy are occupied by 
six electrons, which includes three dangling bonds on the 
carbon atoms and two bonding bonds on the nitrogen 
atom [l| . The nitrogen atom breaks down the Tg sy m- 
metry and gives rise to the trigonal symmetry [bllllof. A 
single NV center have a long-lived spin triplet in its elec- 
tronic ground state whose levels function as a qubit and 
are tunable with an applied magnetic field Cur- 
rently, the applications rely on the spin in the ground 
state, for example, few-qubit networks for simple algo- 
rithms and quantum memories are created by coherently 
coupling NV-center spins to nearby electronic [13, EH 
and nuclear spins [t| EH ; the spin triplet ground state 
also acts as a sensitive magnetic probe of the local envi- 
ronment [20l - [22| . Optical techniques constitute powerful 
tools in quan tum physics and quantum information sci- 
ence [23l |24|. It is known for a considerable time that 
the NV in its singly-charged state is good single-photon 
source [25|. Nowadays, its optical Rabi oscillations have 
been observed 26]. Coherent population trapping us- 
ing optical laser fields has been found in single NV cen- 
ters [27|. And the crucial element for solid-state realiza- 
tion of quantum optical networks has been demonstrated 
by the entanglement between the polarization of a single 
optical photon and a solid-state qubit associated with the 
single electronic spin of a NV center in diamond [28| . 

An important advantage of optical techniques in com- 
parison with microwave manipulation is spatial resolu- 
tion which can selectively address a single NV defect in 
diamond as well as individual qubits. Here, we propose 
an scheme to achieve an qubit rotation by the coherent 
interaction of light and single NV centers in diamond. 
The qubit is encoded in the degenerate Zeeman sublevels 
in the absence of external strain and electric or magnetic 
fields. Due to the spin-orbit and spin-spin interactions, 
an excited-double four-level system can be established in 
the spin-triplet ground state and an orbital-doublet, spin- 
triplet excited state of the NV center, which allow us to 
apply two pump and two Stokes laser pulses. The rota- 
tion of qubit based on the non-Abelian geometric phase 
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is realized by tuning adiabatically the Rabi frequencies 
and phases of the pump and Stokes pulses. 

The paper is organized as follows. In Sec. lU the en- 
ergy levels of a NV center is briefly reviewed and the 
excited-double four-level system is modeled for a NV cen- 
ter with an electronic spin triplet in the ground state and 
an orbital-doublet, spin-triplet excited state. In Sec. IIII1 
we present our method for arbitrary qubit rotation. Sum- 
maries are made in Sec. IIVI 



II. FORMULATION OF AN 
EXCITED-DOUBLE FOUR-LEVEL SYSTEM 

NV centers consist of a substitutional nitrogen atom 
and an adjacent vacancy, which are a naturally occur- 
ring defect in diamond with highly localized electronic 
bound states. The bound states of NV centers are mul- 
tiparticle states of six electrons: five contributed by the 
nearest-neighbor nitrogen and carbon atoms to the va- 
cancy and an additional captured from the bulk. To find 
out the electronic states of NV centers, single-electron 
molecular orbitals are employed, which are built up from 
linear combinations of the dangling sp3 orbitals of one 
nitrogen and three carbon atoms around the vacancy. 
Applying group theoretical arguments to the NV center 
leads to four molecular orbitals {u, v, e x , e^} [l3 - IT5j . Or- 
bitals u and v are totally symmetric but non-degenerate 
which transform according to the one-dimensional irre- 
ducible representation of C% v symmetry, and the two 
degenerate orbitals {e x ,e y } transform according to the 
two-dimensional irreducible representation of C^ v point 
group. The ordering of the orbitals are obtained by 
symmetry and charge considerations of the electron-ion 
interaction, i.e. u orbital has the lowest energy, v is 
the next lowest orbital and doubly degenerate e XtV or- 
bitals are the highest energy. Filling four orbitals with 
six electrons according to the Pauli exclusion principle 
indicates that the ground state consists of four elec- 
trons completely occupying the totally symmetric or- 
bitals and two electrons occupying the remaining or- 
bitals. Although the spin-obit wave functions for the 
ground-state configuration give 3 A 2 , 1 A\ and 1 E states 
due to the antisymmetry property of the total wavefunc- 
tion for fermionic particles, Hund's rules predict that 
the ground state is attributed to the spin-triplet 3 A 2 
state. The excited state configuration can be learned 
from one electron being promoted from the v orbital to 
the e XiV orbitals in the single particle picture 12]. The 
excited-state configuration has an orbital-doublet, spin- 
singlet 1 E state and an orbital-doublet, spin-triplet 3 E 
state with manifold {E x , E y , E' x , E' , A\, A%\ according 
to the irreducible representation of C^y symmetry. The 
degeneracy between triplets and singlets in the ground- 
state and excited-state configurations are lifted by the 
electron-electron Coulomb interaction. Nowadays, it is 
known that the optical transition is associated with the 
3 A 2 ground state and 3 E excited state. The spin-spin in- 
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FIG. 1: (Color online) The energy diagram and selection rule 
of the excited-doublet four-level system associated with the 
ground and excited states in the negatively charged nitrogen- 
vacancy center. 

teraction causes the 2.87GHz zero field splitting between 
the magnetic sublevels mg — and mg = ±1 states of 
the 3 ' A 2 ground triplet state. The spin-orbit interaction 
splits the 3 E excited state into three doublet degenerate 
states, noted by E,E',Ai and A 2 . And the degeneracy 
of the A\ and A 2 states is lifted by the spin-spin inter- 
action. 

The cr + (er_ )-polarized optical pulse drives the transi- 
tion between the ground m s = —1 (+1) state with the Ai 
and A 2 states. Neglecting the <r + (c_)-polarizcd optical 
coupling between the ground m s = 1 (—1) state and the 
E x and E y states and the x (y)-polarized optical coupling 
between the ground m s — state and the E' y (E' x ), a 
excited-doublet four-level system is therefore configured 
in the absence of external strain and electric or mag- 
netic fields, whose energy diagram and selection rule are 
schematically illustrated in Fig. [TJ The ground states 
with m s = ±1 are denoted by |1) and |2), and the ex- 
cited states are denoted by |3) and |4). The transitions 
from states |1) with energy uj\ and |2) with energy u 2 to 
the upper two states with energy L03 and W4 are driven by 
the pump fields and the Stokes laser fields, respectively. 
The pump field is a er_ -polarized bichromatic wave with 
frequencies v p \ and v p2 . And the Stokes field is a cr + - 
polarizcd bichromatic wave with frequencies v s \ and v s2 . 
The Hamiltonian for this system is then 

H = Ul |1) <1| + w 2 |2) (2| + w 3 |3) (3| + uJi |4) (4| (1) 
+ l - (Sl vl e- iv ^ - n p2 e- iv ^) (|1) (3| - |1) (4|) + h.c. 

- % - (firie-*"' 1 * + n s 2e-^ 2t ) (|2) (3| + |2) (4|) + h.c. 

where fl p i and tt S i are time-dependent Rabi frequencies 
where a n phase difference has been introduced in the 
pump field. By assuming that frequencies v p i = v s i, we 
can define two-photon detunings Ai = — w\ — v v \ and 
A2 = v p2 + — UJ4 due to the degeneracy of the ground 
state. In the rotating frame with respect to 
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H = \i) (i\ + («3 ~ Ai) |3) (3| + (w 4 + A 2 ) |4) (4| , 

(2) 
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FIG. 2: (Color online) The norm of probability amplitudes for 
the 3 A2 ground and 3 E excited states as function of time when 
NV center is initially in the ms = — 1 state of the 3 A2 ground 
triplet state, (a) for ms = — 1 ground (blue) and A2 excited 
(red) states, (b) ms = 1 ground (blue) and Ai excited (red) 
states. The norm of the amplitude for ms = ground state 
is plotted in (c) with blue. The superpositions of states E' x , 
E x , E y and E' y give rise to other four states whose norms are 
shown in (c) with red color and (d). All Rabi frequencies are 
taken as 14MHz. And detunings Ai = A 2 = 10MHz. Other 
spin-orbit and spin-spin parameters are taken from Ref. [l4| . 
The time axis is in units of ns. 



the Hamiltonian in the interaction picture is obtained as 
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by employing the rotating-wave approximation. 

To justify that it is valid to neglect optical coupling be- 
tween the ground m s — 1 (— 1) state and the E x and E y 
states, as well as the coupling between the ground m s = 
state and the E' y (E' x ), we take the spin-orbit and spin- 
spin parameters from Ref. [Hj]. Besides, we set all Rabi 
frequencies f2 P j and fl S i have the same value 14MHz, 
and detunings Ai = A2 = 10MHz. Then we numerically 
solve the time-dependent Schrodinger equation for the 
spin-triplet 3 A 2 state and orbital-doublet, spin-triplet 3 E 
state with the system initially in the ms = — 1 state of 
the 3 A 2 ground triplet state. Figure [5] numerically illus- 
trates the norm of probability amplitudes as a function 
of time for (a) ms = — 1 ground (blue) and A 2 excited 
(red) states, (b) ms = 1 ground (blue) and A\ excited 
(red) states. The norm of the amplitude for ms = 
ground state is plotted in Fig.(c) with blue. The super- 
positions of states E' x , E x , E y and E' give rise to other 
four states whose norms are shown in Fig. [U[c) with red 
and (d). It can be found that the probabilities on the 
ground ms — ±1 states and Ai,A 2 excited states are 
much larger than those of other states. Therefore, we can 
eliminate the contribution of the other states and model 
the NV centers as excited-doublet four-level systems. 



III. QUBIT ROTATION 

We now find the eigenstates of Hamiltonian ([3]) . The 
dark states would be desirable for gate operation in or- 
der to remove the dynamic phase shift. Dark states are 
zero-eigenvalue eigenstates of a considered system. The 
characteristic equation of this excited-doublet four-level 
system shows that two eigenvalues of Eq. ^ will be zero 
if we let frequency detunings Ai = A 2 = u>o, Sl p i = il p 
and fl s i — tt s . And these two corresponding dark states 
\Di) and \D 2 ) are coherent superpositions of the bare 
states and have the following expression 
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Here, 
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is orthogonal to the dark state \D\). 

Let us first consider that the applied pump and Stokes 
fields are in resonance with its corresponding atomic 
transition, and two-photon resonance is still maintained. 
In this case, we have wo = 0. The dark state I-D2) is co- 
herent superpositions of the bare states |3) and |4), which 
leads to {Di \ d t \Dj) =0 for i ^ j, i.e. the non-adiabatic 
transition between the two degenerate dark states van- 
ishes. Then, one can completely transfer the population 
from |1) to 1 2) vice versa by the techniques of Stimulated 
Raman adiabatic passage since the pump and Stokes field 
are involved in the dark state \D{). 

When the applied fields are off resonance, i.e., wo 0> 
there is a non-adiabatic transition between these two 
dark states, which plays the important role in single- 
qubit gate operation. Notice that states |1) and \2) span 
the same eigenspace as the dark state |Di) and bright 
state l-Bi), and the qubit described by a|l) + /3|2) does 
not change with time when fl s = Q p = 0, we can project 
the qubit on the dark state \D\) and bright state \B\). 
The projection requires the Rabi frequencies fl p and O s 
to have essentially the same envelopes 



n„(t) =fl(t) cosxe^e^ (t) , 



i<j> iif>(t) 



(6a) 
(6b) 



Here, x an d 4> arc time-independent angles, which specify 
the Euler angle of a vector 

|D) = sinxe"^ |1) + cosxe i0 |2) 

The first dark state and bright state can be written as 

ID!) = e -*iW |D) , = e^W |_ D ) . (7) 
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We denote the Bloch sphere containing vectors |D) and 
|— D) as BS1. The fixed angles \ an d 4> keep vectors |D) 
and |— D) as constant vectors in BS1 during the time evo- 
lution. Hence, any initial state of the qubit can be writ- 
ten as a superposition of |D) and |— D). The mutually 
orthogonal state vectors |— D) and |0) = (|3) — |4))/v / 2 
construct the north and south pole of another unit 2- 
sphere denoted as BS2. The second dark state I-D2) is a 
vector in BS2, whose adiabatic evolution can be written 
as 



\D 2 ) = isin<p (t) |0) + e^W cos ip (t) |-D 



(8) 



where tan ip — 51/ (v2cjo). Since the two photon detun- 
ing is always on, the adiabatic pulses force the system to 
remain in dark states. However, vector I-D2) moves slowly 
with time due to its time-dependent azimuthal and pole 
angles. When the pump and Stokes pulses are adiabat- 
ically turned on, the eigenstate I-D2) evolves smoothly 
from the initial state |— D) to a intermediate state which 
is a superposition state of |0) and |— D). The inverse pro- 
cess, i.e. adiabatically turning of the pump and Stokes 
pulses, brings the system back to its initial state |— D). 
When angles ip and if) make a cyclic evolution, the path 
of the I-D2) vector encloses a patch on the BS2 with a 
solid angle j c , which is the Berry phase. This cyclic evo- 
lution gives rise to a phase change on the state |— D). 
The unitary operator 



U c = e llc |-D) <-D| 



ID) (D| 



(9) 



depicts the effect of the cyclic evolution on the vectors in 
BS1. In the bare basis |1) and |2), the unitary operator 



cos — 1 + 1 sin — n ■ a 
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(10) 



rotates the qubit encoded in states |1) and |2) about the 
axe 

h = (— sin2x cos 2<^), — sin2xsin20, cos 2x) (11) 

by an angle "f c apart from a global phase 7 c /2, where 
a = {(J x ,a y ,a z ) are the pauli's spin operators. Hence, 
arbitrary rotation of the qubit can be accomplished by 
an adiabatic process. 

We now design a rotation of the qubit about x with an 
angle 7r/2 to support the above analysis. The rotation 
about axe a; is obtained by choosing fixed angles \ = 
— 7r/4 and <j> — in Eq.©. For an adiabatic evolution, 
f2 (t) and tp (t) are required to be smooth functions which 
vanish at infinity. We take the envelope and phase of 
Rabi frequencies as follow 
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FIG. 3: (Color online) The norm of amplitudes for state |1) 
(blue dot-dashed line), |2) (red dashed line), |3) (yellow dot- 
ted line) and |4) (green solid line) as functions of time in the 
double-excited model described by Hi in Eq.Q. The hori- 
zontal axe is in units of ns. 



so that the vector \D 2 ) initial at north point undergoes 
an adiabatic evolution along a closed path with a solid 
angle j c = ir/2 in the BS2, where 



77 = — [1 + tanh (at)} 



(13) 



Parameter a characterizes the bandwidth of pulses. Hav- 
ing the designed pulses (|12[) . now we can test our main 
theoretical results (fTU)) with the aid of numerical calcula- 
tions. In order to do the calculation, we first assume that 
the system is initially populated in state |1). Then we 
numerical solve the time-dependent Schrodinger equation 
described by Hamiltonian in Eq. (j3j with loq — 20MHz, 
a = 1MHz. In Fig. [3J we plot the norms of ampli- 
tudes for all four states as functions of time. It can be 
found that there are no population in the two excited 
states (yellow dotted and green solid lines) before and 
after applying the pulses, however the amplitude of the 
initial state |1) (blue dot-dashed line) has changed from 
1 to 1/V2, and the amplitude of the degenerate state 
1 2) increases from to 1/V2. Therefore, a rotation of 
the qubit about x with an angle 7r/2 has been obtained. 
The non- vanishing amplitude of the double excited states 
during the adiabatic evolution depicts the non-adiabatic 
transition between two dark states. 

Notice that the optical transitions in NV centers are 
original associated to an orbital-doublet, spin-triplet ex- 
cited state and the spin-triplet ground state, the actual 
operation U is presented by 9 x 9 matrices. To check 
how well our theory works for actualy nine-level sys- 
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terns, we use the operator fidelity [2 

as a measure of how accurate the ideal gate Uc can be 
achieved. Taking the spin-orbit and spin-spin parame- 
tersf from Ref.pjj], and letting the Rabi frequency am- 
plitudes, the pulse widths and the detuning are same as 
those in the above, we numerical solve the equation for 
the unitary evolution operator with the initial condition 
U (to, to) = I- The actual operation U is obtained at the 
time t f that the applied pulses have turned off. Then we 
obtain the operator fidelity on the order of 99.9%. 
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IV. CONCLUSION 

In conclusion, we have shown that it is possible to es- 
tablish an excited-double four-level system in the spin- 
triplet gound state and an orbital-doublet, spin-triplet 
excited state of the NV center due to the spin-orbit 
and spin-spin interactions, which allow us to apply two 
pump and two Stokes laser pulses. It is found that two 
dark states exist when both the pump and stokes fields 
maintain two-photon resonance and the intensities of two 
pump (Stokes) fields take the same values. Consequently, 
we propose a rotation procedure by encoding the qubit 
in the ground degenerate states. We first map the in- 
formation of the qubit on one dark state and its corre- 
sponding bright state defined by the pump and Stokes 
pulses. Then we adiabatically tune the Rabi frequencies 
and phases of the pump and Stokes pulses by keeping 
nonvanishing two-photon resonance so that the NV cen- 
ter undergoes a cyclic evolution. The adiabatical process 



first transfers some components of the bright state to 
the double excited states, afterwords, bring these com- 
ponents back to the bright state accompanied by a phase 
shift 7c. Therefore, a rotation of the qubit is acheived 
due to the non-Abelian geometric phase produced by the 
nonadiabtic transition between the two degenerate dark 
states. The axis and angle of rotation is determined by 
the parameters of the laser fields. 
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